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Abstract 

The temperatures in different zones in tlie world do not show significant changes due to El- 
Nifio except when measured in a restricted area in the Pacific Ocean. We find, in contrast, that 
the dynamics of a climate network based on the same temperature records in various geographical 
zones in the world is significantly influenced by El- Nino. During El-Niho many links of the network 
are broken, and the number of surviving links comprises a specific and sensitive measure for El- 
Niiio events. While during non El-Niho periods these links which represent correlations between 
temperatures in different sites are more stable, fast fluctuations of the correlations observed during 
El-Niho periods cause the links to break. 
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Networks are often used to formulate the dynamics of complex systems that are built 
from many interacting components (see e.g. PQ[2])- While in some cases the representation 
of the system as a network is obvious and the nodes and links are identified directly (e.g. 
cables connecting computers in a computer network) [3], there are cases in which the process 
that couples the individual interacting components is more complex and a link is guessed 
by tracking similarities in the dynamical behavior of two nodes [1]. 

Even when the usual traces of dynamics of interacting nodes on a network, such as 
partial synchronization, clusters with correlated dynamics, oscillatory synchronization [S], 
and phase slips [B], are evident, the mission of designing a generic tool that reliably extracts 
information about the network structure from measurements of the dynamics of nodes is 
still far from being accomplished. 

In this Letter we develop a method for generating climate networks, which is suitable 
for tracking structural changes in these networks (dynamics of a network). These changes 
correspond, in our case, to strong climate changes due to El-Nino. We find that networks 
constructed from temperature measurements on different sites in the world are changed 
dramatically during El-Niiio events in a similar way. These structural changes are seen even 
for geographical zones where the mean temperature is not affected by El-Nino. 

We analyze daily temperature records taken from a grid (available at [7|) in various 
geographical zones (shown in Fig.[T|. To avoid the trivial effect of seasonal trends we subtract 
from each day's temperature, the yearly mean temperature of that day. Specifically, if we 
take the temperature signal of a given site in the grid to be T^{d), where y is the year and d 
is the day (ranging from 1 to 365), the new filtered signal will be T^^d) = T^i^d) — ^ J2y T^^d) 
(where N is the number of years available in the record) [8]. 

We compute for a time shift r G [—Tmax, Tmax] days for each pair of sites / and r on the grid, 
their cross correlation function Xi,^,{t > 0) = {Ti{d)Ty{d + r))^ and Xf.^{T < 0) = Xr,ii'^ > 
0). The correlation strength of the link is chosen to be Wj'^ = MAX{Xf^)/STD{X^^), 
where MAX and STD are the maximal value and the standard deviation of the absolute 
value of Xi^ in the range of r, respectively [9]. The time shift at which Xf,, is maximal 
is defined as the time delay. Up to here, the prescription is similar to other methods (see 

e.g. m )■ 

From reasons that will become clear later, we are able to set a physical threshold Q so that 
only pairs /,r that satisfy Wi^. > Q, are regarded as significantly linked. Mathematically 
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this can be represented by the Heaviside function 0( follows, 



Plr = QiWyr-Q)- (1) 

Some of the elements of the matrix p may blink as a function of y, i.e., appear and 
disappear. Thus, even though there are many pairs l,r that their correlation values W^f^ in 
a specific year are above Q, some of these Wl'^ are sensitive to the choice of the beginning of 
the period y, and to noise . In the present work we choose to discard the question of which 
pairs comprise the static network, and concentrate in the structural changes of the network 
over time. Blinking links seem to be a signature of structural changes, so we distinguish 
between them and the more robust links, that are stable during larger time periods. 

In the next step we examine if a currently existing link pf ^ existed in the earlier periods of 
the network. To accomplish this we define a new matrix, which takes into account previous 
states of the links in the last k states of the network. We define a new matrix Mf^ which 
counts the number of times a link appeared before continuously (without a blink): 

K = T. n pZ- (2) 

n=0 m=y—n 

A link {l,r) in the current network p appeared k times in a row before (including its 
current appearance) iff M/^^ > k. These links represent long lasting relations between 
temperature fluctuations in the zone. Counting them enables us to distinguish between the 
two qualitatively different groups of links, blinking links which are removed, and robust links 
which we include in the network. The number of links that exist in our network depends on 

N N 

nk{y) = Y: E 0{Mlr-k + l). (3) 

1=0 r=l+l 

Where k is the number of times a link has to survive in order to be included in the network 
(according to Eq. (|2])), and is the total number of links in the geographical zones. The 
summand in the rhs of Eq. ([s]) represents the network matrix. In the current work we chose 
the y resolution (the jumps between two subsequent dates represented by y) to be 50 days, 
k = 5 and the threshold Q = 2. 
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FIG. 1: The four geographical zones used for building the four climate networks studied. The dots 
represent the nodes of the network. The rectangular geographical zone inside zone 1 shows the 
standard basin for which El-Niiio effects on temperature and pressure is significantly observed (see 
Fig.|2}i) 

The results shown below are not sensitive to the choice of k. However, choosing too 
large k values reduces the number of surviving links significantly, and therefore eliminates 
much of the effect. Choosing too small values of k, on the other hand, does not enable the 
elimination of blinking links, and therefore causes Uk (y) to be more noisy, but the significant 
effect of breaking links is still evident. 

We chose four representative zones around the globe, as shown in Fig. [1} Our networks 
are built from measurements of temperatures close to sea level (networks Aj), and from 
measurements on a 500mb pressure level (networks Bj), on a grid of 7.5° resolution. The 
measurements are taken for the years 1979-2006, for which 8 known El-Nifio events have 
occurred . 

In Fig. [2^ we show the effect of El-Nifio as 8 main extreme values on the two standard El- 
Niiio indices (based on temperature and pressure measurements), measured in the standard 
basin region inside zone 1 of Fig. [l] [13]. In the top of each of the eight panels in Fig. |2}d,c we 
show the mean temperature anomaly over the whole corresponding zone defined in Fig. [1} 
It is clearly seen that compared to Fig. [2^ the El-Niiio effect on mean temperature in all 
zones becomes very weak and almost cannot be detected except at zones Ai and Bi (which 



4 



(a) 





0.5 

-0.5 



(b) 



4250 
4200 
4150 - 




19% 1^53 — Tm 1953 — 35D — — 1*3 — m>_ 1553 — W) W 



1 

0.5 

-0.5 
-I 

2]000 

20000 
19000 
18000 
17000 
16000 



4280 
4270 
4260 
4250 
4240 



r- 



790( 
7S0f 
7701 
760( 



.1500 
.3450 
3400 



V7 



4 



ii 



1980 1985 1990 1995 2000 2005 1980 1985 1990 



FIG. 2: (a) Mean sea surface temperature (left) in the standard basin shown in the rectangle inside 
Ai in Fig. [l] (NIN03), and the difference in sea level pressure (right) between Tahiti and Darwin 
(SOI), both are standard indices for El-Niiio (see e.g. [H]). (b) and (c) (color online) The upper 
curves represent the temperature anomaly series in zones (b) Ai, A2, A^, A4^ and (c) -61,-62,-63,-64. 
The lower curves present (y), namely, the number of links that survives in the network as a 
function of time, for these same zones. In zone -64, the graph for (y) is completely flat, when 
Q = 2. For this zone we show the cases of Q = 2.4 (middle curve) and Q = 2.5 (lower curve) [12] 
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include the standard basin). In most zones, El-Nino effect on the temperatures are masked 
by noise of the same order. In mark contrast, when measuring the number of links {y) 
in the climate network we observe (in the bottom of Figs. |2|d,c) a significant effect of El- 
Nifio in all four zones. This is represented by the sharp fall of rik (y) in most of the times 
of occurrence of El-Nino. In zones Ai,A2,A^ between 5 and 6 out of the 8 events can be 
observed while in zone (in the north) only the strongest El-Nino is seen clearly. Because 
of the close occurrence of the 3 events in the early 90 's it seems that two of the minima 
of rik (y) overlap. It is also notable that in the zone surrounding the El-Nino basin (zone 
Ai) the minimum due to the 2002 event cannot be inferred. However in the Bi (Fig. ^) 
network's (y) series the can be seen. In zone A^ it is seen that even when the influence 
of El-Nino on the mean temperature in the zone is not seen at all, the largest El-Niiio event 
in 1997 as well as the event of 2005 can clearly be identified by a sharp minimum in Uk [y). 

In zone Bi (Fig. ^) the main events are also clearly seen, but the 3 events in the early 
90's overlap completely. In network B2 all events, except the 2002 event, can be observed. 
In network B^ only the 1997 event is observed reliably. In network i?4 the traces of El-Nifio 
can also be speculated if one sets the threshold to a higher value. 

The choice of the threshold Q = 2 is not arbitrary. When observing the probability 
density function of W (Fig. [3]) it is clear that for non El-Nino time regimes, W = 2 is 
actually the minimal value that exists. It therefore appears that choosing this threshold 
makes the network very sensitive to El-Nino events while remaining insensitive to other 
changes in climate. The reason is that the distribution of Wl'^ tends to typical lower values 
of VFfj, during El-Niiio, as can be clearly seen in Fig. |3]^A-E). A remarkable property of this 
softening is that the lower limit of the distribution drops from being close to 2 to some 
significantly lower value. Changes in climate around the world due to El-Nino events thus 
share a unified property of the correlation pattern, which can be tracked in a reliable way 
by the number of surviving links n/c (y) in the climate network. 

How susceptible are the links that will break during El-Nifio? As can be seen in Fig. 
[3]^F-J), the links in the fully connected network that will later break have distribution of 
correlation strengths very similar to that of the links that will survive. Only a slight tendency 
of weaker links to break can be observed. However, when looking at the distribution of the 
time delays (Fig. |4]), there is a clear difference between links that survive and links that 
collapse. Links that collapse have a much broader distribution of time delays, and therefore 
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FIG. 3: (color online) A-E: The distribution of W^'^. in zone 1 (see Fig.|l|. The five rows correspond 
to the five El-Nifio events pointed with an arrow in Fig.|2^, which are the most noticeable on zone 1. 
The curves with a peak to the right (blue) describe the distribution in a reference time, where the 
influence of El-Niho on the network is not seen. The curves with a peak to the left (green) describe 
the distribution for five El-Niho events. F-J: Two distributions of Wj^^ in the reference point (blue 
curves of graphs A-E). The curves slightly shifted to the left (black) are the distributions of the 
correlation values of the links that later on will disappear due to one of the El-Niho events , and the 
curves slightly shifted to the right (red) are the distribution of the ones that will remain connected 
in the considered event. 
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FIG. 4: (color online) Distribution of time delays. A-E, Lower (green) curves : The distribution 
of time delays of all links in the reference time, with no signs of El-Niho. Upper (blue) curves: 
Time delays distribution during five El-Niho episodes. F-J : Distribution of time delays of two 
groups of links during the reference time. Upper (red) curves - distribution of time delays of links 
that remain during El-Niho. Lower (black) curves- distribution of time delays of links that break 
during El-Niho. 
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have usually delays of more than a few days. Links with shorter time delays have higher 
probability to survive. 

An alternative approach that yields a similar effect can be obtained by counting the 
number of elements of pf^ that during the periods {y— 2, y—l, y} (two subsequent dates y,y+l 
are separated by 50 days) posses for pf^. the sequence {1, 0, 1} or the sequence {0, 1, 0}. Such 
a treatment which tests directly the relation between El-Nino and blinking of links yields 
results very similar to that given by (y). The blinking correlations events during El-Nifio 
may arise from a gradual breaking of indirect paths between two nodes / and r, or from a 
sensitivity to the choice of Q. On either case these events are a direct and unified mark of 
an aggressive structural change, that until now was described as various climate responses 
in various geographical zones. A more resolute theoretical problem of reconstructing an 
accurate network structure from experimental observations of its nodes' dynamics (in the 
current paper - the temperature ) remains largely open (see related discussions in [HI [15], 
)• 

In summary, we have developed a method which enables one to follow large changes over 
time of an underlying network structure by observations of fluctuations in the correlations 
between nodes. The method tracks blinking links that appear and disappear in a short 
time, and assumes this behavior to be due to structural changes. Tracking the changes in 
the network of temperatures in several zones in the world reveals a deep violent response to 
El-Nino even in zones and heights where the mean temperature level is not affected. The 
links that break during El-Niiio are mostly links that have large time delays. 
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